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ABSTRACT 

Lukens , Larry Allan, Ph.D., Purdue University, 
March 1971. An Experimental Investigation of Electric 
Field Intensity and Wall Heat Transfer for the Heating 
Region of a Constricted Arc Plasma. Major Professor: 
Frank P. Incropera. 

Electric field intensity and total wall heat 
transfer measurements have been performed for the inlet, 
asymptotic, and anode regions of a laminar, atmospheric 
cascade arc. In addition to the diagnostics required 
for these measurements, a unique radiation gage has 
been developed and used to obtain radiative wall heat 
flux data. While the majority of the measurements have 
been made for an argon arc, some data are reported for 
helium and nitrogen arcs . 

Comparisons between the experimental data and 
available equilibrium and nonequilibrium numerical 
solutions are presented in an attempt to ascertain the 
adequacy of single-fluid models to accurately describe 
overall arc energetics. For the most part, the equili- 
brium theory predictions for the asymptotic region in an 
argon arc agree with the experimental results to within 
the combined experimental and theoretical uncertainties. 
However, there is a small, but discernable nonequilibrium 
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effect at both the lower and higher arc currents in argon. 
In contrast, the nonequilibrium theory predictions are 
in excellent agreement with the experimental results over 
the entire range of asymptotic conditions investigated in 
argon. Further comparisons are presented which indicate 
that atmospheric helium arcs are characterized by a signi- 
ficant departure from local thermodynamic equilibrium. 

On the basis of comparisons between experiment and 
theory for the anode region, it is shown that anode heat 
balance models are ill-suited for providing design es- 
timates for the total heat flux to the anode of the arc 
constrictor . 
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CHAPTER I INTRODUCTION 

1.1 General Objectives of This Study 

Owing to interest in space flight technology, MHD 
systems, and thermonuclear energy conversion, considerable 
effort has been devoted to the study of high temperature 
gas flows. Of particular interest to this study is that 
portion of this effort which deals with a device which has 
long served as a versatile tool in generating high tempera- 
ture gases for both industrial and laboratory use - the 
plasma arc generator. 

Designed initially by Maecker^^, the operation of 
the Gerdien, or constricted tube, arc generator is based 
upon the conversion of electrical to thermal energy. In 
the most common configuration (Figure 1.1), this conversion 
is accomplished by passing a gas through an electric arc 
which is established by applying a d-C potential between 
the electrodes. 

The constricted tube arc produces a continuous and 
well-defined plasma flow which, in addition to its use in 
certain technological applications, is particularly useful 
for the experimental determination of high temperature 
gas transport properties and the development of high 
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temperature diagnostic techniques. Irrespective of its 
use, however, the constricted tube arc is distinguished 
by extremely high heat fluxes to the containing walls and 
these walls are often destroyed if sufficient cooling is 
not provided. Accordingly, an understanding of the elec- 
trical and thermal characteristics of arc constrictors 
is important to the designer as well as to those with a 
general interest in arc physics. The objective of this 
study is to enhance this understanding through an experi- 
mental investigation of the electrical and thermal effects 
in an atmospheric cascade arc and through an interpretation 
of these effects in terms of available theoretical studies. 
While the major portion of the experimental work has been 
performed for an argon plasma, some measurements have 
also been made for nitrogen and helium arcs. 

The regions of interest for the cascade arc are 
depicted schematically in Figure 1.1. The arc column, 
which is a composite of individually water-cooled and 
electrically isolated copper discs, is divided into three 
distinct regions: inlet, asymptotic, and anode. The gas 
flow in each region is subjected to an externally applied 
electric field, and the inlet and asymptotic regions are 
generally bounded by non-current-carrying surfaces. In 
contrast, the anode surface is electrically conducting. 
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Figure 1.1 Designation of Arc-Heating Regions 
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The inlet region is further distinguished by a significant 
axial variation of the flow variables, whereas the asymp- 
totic region is characteristic of thermally fully-devel- 
oped flow. As a matter of convention, the axial location 
at which the mean (mass averaged) enthalpy attains 90% 
of its fully-developed value is chosen to mark the be- 
ginning of the asymptotic region. 

1.2 Previous Investigations of 
Plasma Tube Flow 

Although the literature devoted to theoretical study 

of the inlet and asymptotic regions is extensive, few 

experimental studies have been performed for purposes of 

determining the overall arc energetics. Of particular 

scarcity are experimental studies pertaining to the inlet 

[21 

region of the water-cooled duct. Skifstad performed 

electrical potential and wall heat flux measurements in 

the inlet region of a helium arc, and investigators at 

T 3-51 

the NASA-AMES Research Center have made similar 

measurements in air and nitrogen arcs. However, the range 
of experimental conditions (gas flow rate, arc current, 
and duct geometry) investigated in these studies differ 
significantly from those of this investigation. Conse- 
quently, meaningful comparisons cannot be made between the 
data of this study and the experimental results reported 



5 



r 6 7 1 

by these investigators. Similarly, Pfender, et al. ' 
report results for the local wall heat flux in the inlet 
region of an argon cascade arc; however these measurements 
were for turbulent flow and were used primarily for the 
support of anode heat transfer studies. 

A greater number of experimental investigations 
have been performed for the asymptotic region of the 
arc L . While these studies are primarily concerned 

with the use of an inversion method to determine plasma 
transport properties, they necessarily include measure- 
ment of the arc-column electrical characteristics and 

are therefore germane to this study*. In addition, Morris, 
[ 14 ] 

et al. present asymptotic wall heat flux data and 

Emmons and Maecker^^ report on the percentage con- 

tribution of radiation to the total energy balance. 

In view of the excellent reviews given by 
Skif stad and, more recently, by Bower^^, the present 

discussion of theoretical investigations is limited to 
those studies which are to be used for comparison with 
experiment in Chapter IV. 

The approximate analytical solution by Stine and 
[181 

Watson 1 was one of the first reported theoretical 
treatments of the inlet region of the cascade arc. This 



* The inversion method of property determination will be 
discussed further in Chapter IV. 
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solution is based upon a reduced form of the governing 
energy equation, obtained by neglecting radiation, radial 
convection, and viscous dissipation. In addition, the 
Stine-Watson model assumes a uniform mass flux throughout 
the constrictor and uses linearized relationships between 
the enthalpy and related properties (T , k/C , and a). The 
primary weakness of this model is its failure to properly 
account for the variation of the transport properties 
with temperature and, particularly when applied to the 
higher temperature argon and nitrogen arcs, its failure 
to account for radiative losses. 

To remove the need for the above assumptions, de- 
tailed numerical solutions for the inlet region have been 

[ 19 ] 

carried out by Watson and Pegot and by Bower and 

Incropera . Of the two, the latter will be used for 
the comparisons of this study since it gives greater con- 
sideration to pertinent thermal effects. Specifically, 
Bower and Incropera formulate a rigorous equilibrium model 
for the laminar, steady, and axially symmetric (without 
swirl) flow of an atmospheric argon plasma through the 
heating region of a constricted tube arc. An implicit 
finite-difference scheme 1 is then used to solve the 
resulting conservation equations (mass, momentum, and 
energy) consistent with the thermal and caloric equations 
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of state. Ohm's law, and the temperature dependence of the 
required thermodynamic and transport properties. Implicit 
in the final form of the equations are the boundary layer 
approximations and the assumptions of negligible radial 
current density, quasi-neutrality, and local thermodynamic 
equilibrium. The results of this numerical solution are 
presented in the form of wall heat transfer and shear 
stress (as predicted from energy and momentum balances) and 
correlations between the mean Nusselt number and an Ohmic 
heating parameter. 

As a result of the evolution of the aforementioned 

numerical solutions, it would appear that the theoretical 

treatment of the inlet and asymptotic regions of arc 

constrictors has reached a practical limit. However, the 

r 13 1 

variable tube diameter experiments performed by Emmons 

and, more recently, the spectroscopic measurements con- 

(211 (221 
ducted by Bott and by Incropera and Giannaris 

provide strong experimental evidence of the presence of 

nonequilibrium in confined plasmas. Accordingly, the 

validity of the local thermodynamic equilibrium assumption 

which is implicit in the single-fluid model is subject to 

question . * 



A comprehensive discussion of local thermodynamic 
equilibrium in confined plasmas is given in references 

and ' 23 h 



[ 22 ] 
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In an attempt to account for the nonequipartition 

of species translational energies (thermal nonequilibrium) , 

[241 

Pytte and Winsor formulated a nonequilibrium theory 

for the asymptotic region of a hel.ium arc. This attempt, 
however, is at best, semi-empirical, for its solution 
relies upon experimentally determined arc-column charact- 
eristics. In addition, the Pytte-Winsor model suffers 
from use of the Saha (equilibrium) equation to compute 
the plasma composition. 

While accounting for both thermal nonequipartition 

f 25 1 

and chemical nonequilibrium, Clark has formulated a 

rigorous multifluid (electrons, ions, and neutral atoms) 
model for the steady, laminar flow of a nonequilibrium 
argon plasma in the inlet and asymptotic regions of a 
constricted tube arc. After applying the boundary layer 
approximations, the governing set of species and global 
conservation equations, along with the thermally perfect 

gas equation of state and Ohm's law, are solved numerically 

• . . [ 2 5 1 

using finite-difference techniques . Of particular 

interest to this study are the predictions for the asymp- 
totic electric field intensity, total wall heat flux, and 
ratio of radiative to total wall heat flux which result 
from the nonequilibrium solution. 

The need for extended operating times for the most 
severely heated component of the arc constrictor, the 
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anode, has prompted considerable research on the complex 
energy exchange mechanisms occurring at this surface. 
Energy transfer models have been proposed in which the 
various contributions to the total anode heat flux are 
delineated. For the most part, these models are based on 
a heat balance across the anode surface 1 ' ' * . In 

addition, extensive studies of the electrical characteris- 
tics of anodes have been conducted by researchers at 
Dartmouth . For purposes of comparison with the 
anode heat transfer measurements of this study, primary 
emphasis will be placed upon the widely used model which 

r 6 1 

is described by Pfender, et al. 1 1 . 

1.3 Scope and Specific Objectives of This Study 

The primary objective of this study is the measure- 
ment of the axial distribution of the electric field 
intensity and the total wall heat flux for the inlet, 
asymptotic, and anode regions of an atmospheric argon 
cascade arc. These measurements are made for a wide 
range of laminar flow operating conditions. Additional 
measurements are made to determine the ratio of radiative 
to total wall heat flux for the asymptotic region of an 
argon arc**. An important secondary objective includes 

* A review of the subject of anode heat transfer for 
various electrode geometries is given by Eckert and 

Pfender . 

** A review of the literature pertinent to this subject is 
presented in the Appendix. 
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comparison of the asymptotic wall heat transfer and electric 
field intensity measurements with the experimental data of 
other investigators and with the analytical and equili- 
brium and nonequilibrium numerical solutions described 
previously. With respect to the latter, particular atten- 
tion is given to ascertaining both the validity of existing 
Nusselt number correlations and the adequacy of single- 
fluid models for the description of arc energetics. Pur- 
suant to gaining additional information on anode physics 
and the suitability of anode heat transfer correlations, 
further comparisons are made with existing anode heat 
transfer models. The measurements are also of value in 
providing an improved understanding of the development 
length characteristics of the constricted tube, laminar 
arc. In order to ascertain the influence of operating 
gas, additional experimental results are presented for 
nitrogen and helium arcs. 

This study, which is the first serious attempt to 
experimentally verify existing equilibrium and nonequili- 
brium theoretical predictions of thermal parameters for 
the arc heating region, is intended to provide information 
useful to the arc designer and to contribute to an im- 
proved understanding of arc physics and flow modeling. 
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CHAPTER II ARC GENERATOR AND DIAGNOSTIC TECHNIQUES 

II. 1 Arc Generator and Supporting Systems 

A well-defined, stable plasma flow is achieved 

through use of a constricted tube plasma arc generator 

[ 34 1 

manufactured by Creare Inc. (Figure II. 1). The 

constrictor is formed by stacking together a series of 
individually water-cooled copper segments (discs) , each 
of which is 0.632 cm wide with a 1.0 cm central bore. Each 
segment is electrically insulated from adjoining segments 
by cool O-ring seals which further provide for alignment 
and hermetic sealing of the duct. The lengths of the 
arc heating and field-free regions may be varied according 
to the anode location and the number of segments used to 
form the stack (up to a maximum of 40) . The thoriated 
tungsten cathode is water cooled and may be precisely 
centered with respect to the initial nozzle-shaped seg- 
ment through micrometer adjustment. The cathode-nozzle 
geometry provides for a high degree of flow symmetry which, 
along with water cooling, provides for extended cathode 
life. If the tungsten tip of the cathode is periodically 
(every 20 hours) re-shaped on a lathe such that the proper 
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Figure II 
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cathode-nozzle geometry is maintained, a single cathode 
will provide approximately 100 hours of operation. 

The supporting systems which are necessary to sus- 
tain the arc generator are illustrated schematically in 
Figure II. 2. Electrical power is supplied by a D-C 
motor-generator set capable of delivering 75 kw at 300 
amperes. Stable operating conditions are achieved through 
use of two, step-wise variable ballast resistors which 
are connected in series with the power supply and arc 
generator. The electrical system will allow for the use 
of either r-f discharge or prod starting procedures. 

The main coolant system supplies cooling water to 
the plasma generator coolant manifold and to the main 
ballast resistor. The manifold in turn distributes cooling 
water to the cathode, the starting prod, those segments 
which do not serve as calorimeters, and to the additional 
downstream ducting required to cool the gas for discharge. 

A secondary coolant system is employed to reduce flow 
demands on the main coolant system and to provide a 
means of selectively monitoring and controlling the 
coolant flow rates through those duct segments for which 
wall heat transfer measurements are desired. This system 
also supplies cooling water to the auxiliary ballast 
resistor. 
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Figure II 
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Receiving its supply from commerical gas cylinders, 
the gas system facilitates control and measurement of the 
operating gas flow rate. Stringent purity requirements 
are maintained throughout the gas system including use of 
research grade gas, prior ultrasonic cleaning of all 
components exposed to the plasma stream, and continuous 
outgassing of the system prior to operation. 

II .2 Particular Operating Procedures 
and Limitations 

Despite the aforementioned purity requirements , the 
level of duct contamination resulting from arc operation 
has steadily increased over the past two year period. 
During the earlier stages of this research, when the major 

t 

portions of the electric field intensity and total wall 
heat flux data were acquired, there were no visible indi- 
cations of surface contamination, even after extended per- 
iods of operation (up to 40 hours). However, during the 
latter phases of the research, when the radiative wall 
heat flux data were acquired, the level of contamination 
continually increased in spite of special precautions that 
were taken in addition to those listed above. 

The contamination is evidenced by a thin deposit 
which coats the duct walls in a distinct and repetitive 
pattern. The deposit, which appears to be a dark oxide 
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layer, forms on the leading (upstream) edge of each of the 
segments in the arc heating region and diminishes in the 
axial direction, leaving the downstream half of the seg- 
ment surfaces free of visible contamination. The deposit, 
which at times has a bluish cast, has not been observed in 
the field-free region, not even on the segment immediately 
downstream of the anode. 

Similar deposits have been observed by other in- 

r 13 32 331 

vestigators ' ' ; however opinions differ concerning 

the effect of surface contamination on experimental results. 

[ 13 ] 

Emmons claims that the overall arc characteristics are 

virtually unaffected by small amounts of contamination. 

[331 

On the other hand, Runstadler reports a reduction in 

the electric field intensity for helium of some 10% when 
data are acquired under contamination-free conditions. 
Throughout the course of this study, both electric field 
intensity and total wall heat flux measurements have been 
consistently reproducible to within their respective 
experimental uncertainties, irrespective of the level of 
contamination. None the less, it was assumed initially 
that contamination would render the proposed radiation 
gage completely useless, since the dark deposits would 
significantly alter the optical properties of the gage 
surface. While useful data were acquired from the 
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radiation gage, the contamination level within the duct 

still prevented the gage from being used as it was 

originally intended (see the Appendix) . 

[ 32 ] 

As reported by McKee, et al. , the difficulties 
associated with obtaining reproducible measurements of 
anode characteristics are most likely due to the pro- 
nounced affect that surface contamination has on electrode 

performance. In connection with this problem, a number 

[27 32-351 

of investigators L ' 1 have devoted considerable 

effort to the elimination and investigation of surface 

contamination within constricted tube arcs. In summarizing 

[331 

work performed in this area at Dartmouth, Runstadler 
suggests a number of procedures which may be used to re- 
duce the level of contamination within the duct. These 
recommended procedures include use of ultra-pure grade 
gas, elimination of possible cathode contaminants through 
use of a sonic orifice, outgassing of the duct surfaces, 
and nickel plating of all copper surfaces exposed to the 
plasma stream. The plating procedure was found necessary 
in order to eliminate outgassing from the copper surfaces 
during arc operation. 

Similar, yet considerably less involved, measures 
were taken in this study in an attempt to eliminate the 
contamination of the duct surfaces. Ultra-pure grade 
(as opposed to commercial grade) gas was used; however 
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its use did not result in a significant reduction in 
surface contamination. This is not to say, however, that 
use of commercial grade gas will not result in surface 
contamination. After a thorough cleaning of the gas 
system failed to effect any change in the degree of sur- 
face contamination, the inlet portion of the gas system 
was virtually eliminated by connecting the gas cylinder 
directly to the plasma generator (a short piece of plastic 
tubing was connected between the gas inlet to the plasma 
generator and the gas regulator and the regulator was 
connected directly to the gas cylinder). In addition, 
the gas regulator (a two-stage, metal diaphram type) was 
completely disassembled, cleaned ultra-sonically with 
acetone and isopropyl alcohol, and allowed to degas at 

_ g 

10 mm Hg for a period of 24 hrs. Several hours prior 
to starting the plasma generator, the combined generator 
and gas system was evacuated to approximately 40 microns 
and the metal components, particularly the duct segments, 
were heated to accelerate degassing. Unfortunately, none 
of, or any combination of, these procedures led to any 
significant reduction in the level of surface contamina- 
tion. In addition, the possible influence of prod and 
r-f discharge starting techniques was investigated, as 
were the procedures used to mechanically and chemically 
clean the duct segments. Special care was taken in 
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positioning the cathode within the initial nozzle segment 

so as to minimize possible contamination from erosion of 

the tungsten cathode tip. Again, there were no per- 

ceptable changes in the level of surface contamination with 

any of these additional measures. 

On the basis of these results, it is tentatively 

concluded that the surface contamination is an oxide 

layer which is formed by a chemical reaction between the 

copper duct walls and the gases which escape from these 

walls during operation of the arc. This conclusion is 

r 3 3 ] 

supported by the observations of Runstadler concerning 

the nickel plating of the copper surfaces, and also by 
the experiments performed in this study with the thin 
film aluminum surface on the radiation gage (Appendix) . 
Further support of this conclusion is provided by the 

[ 32 ] 

x-ray diffraction analysis conducted by McKee, et al. 
which identified the surface contaminants as copper-oxide 
and carbon. While not fully understood, the form of the 
axial distribution of the surface contamination on seg- 
ments in the arc heating region and the absence of surface 
contamination in the field-free region indicate that 
certain electrical effects must also influence the 



contamination . 
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In addition to experiencing contamination problems, 
some difficulties arose in connection with generating He 
and N 2 arc plasmas. The maximum current-carrying capa- 
bilities of the starting prod (75 amp) and r-f discharge 
(50 amp) circuits, while adequate for initiating argon 
arcs, are insufficient to permit starting in either 
nitrogen or helium. Accordingly, a special pre-start 
procedure was devised. This procedure involves use of a 
dual gas supply manifold whereby both Ar and either N 2 or 
He can be supplied to the arc generator simultaneously. 
Following initiation of a stable, 50 amp Ar arc, the 
proportion of N 2 or He to Ar is slowly increased from 0 to 
100% while the open-circuit voltage is steadily increased 
to maintain a stable operating condition. While this 
procedure is adequate for initiating N 2 and He arcs, the 
maximum obtainable open-circuit voltage of the power 
supply is nonetheless insufficient to sustain stable 
operation of these gases over as wide a range as is 
possible with Ar. 

The restrictions which are placed upon both high 
and low current operation of the arc generator are a re- 
sult of certain physical limitations of the power supply 
and external ballast resistor. Specifically, high current 
operation is limited by the maximum open-circuit voltage 
capability of the power supply. In contrast, low current 
operation is restricted by the lower limit of electrical 
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stability which, in turn, is a function of the available 
ballast resistance as well as the maximum obtainable open- 
circuit voltage of the power supply. 

These limits of operation can be explained by con- 
sidering a set of hypothetical (yet typical) load line 
characteristics for an arc generator and its associated 
D-C power supply (Figure II. 3). As indicated by the 
electrical circuit which accompanies Figure II. 3, the 
power supply load line characteristic is a plot of the 
voltage available to the arc generator versus the available 
arc current. Generally, the power supply characteristic is 
represented by a straight line which intersects the vol- 
tage axis at the value of the open-circuit voltage and 
which has a slope equal to the negative of the sum of 
the constant internal resistance of the power supply (R s ) 

and the variable external ballast resistance (R, ) . The 

b 

arc generator load line, on the other hand, represents a 
plot of the current conducted by the arc versus the 
voltage required to sustain that current. On the whole, 
the arc generator is characterized by a wide spectrum of 
load lines which result from various combinations of the 
operating gas, the operating gas flow rate (m) , the 
operating pressure (p) , the overall (cathode to anode) duct 
length (L) , and the duct radius (R) . 



Increasing 
m and/or L 
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Figure II. 3 Hypothetical Load and Supply Characteristics 
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The actual operating state (for a specified gas, 
ih, p, L, and I) of a constant diameter arc generator is 
uniquely determined by the point of intersection of the 
load (arc generator) and supply characteristics. The 
degree of electrical stability which exists for a given 
operating state is a function of the relative slopes of 
the two load lines at their point of intersection. The 
criterion that must be satisfied in order to achieve 
stable operation is given by 



(if) " (fl) - 0 

LOAD SUPPLY 



(II. 1) 



which, in the case of a linear supply characteristic, 
reduces to 



(£E) 

*dl ; 



+ R 



LOAD 



+ R b ' 0 



(II. 2) 



It is noted at this point that, while Equation II . 2 
establishes, in principle, the criterion for stable opera- 
tion, in practice the slope of the load characteristic must 
always be greater than the slope of the supply character- 
istic if stable operation is to be realized. That is. 
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LOAD 



+ R + R, >> 0. 



s 



(II. 3) 



Since R is constant and since 

s dl 



LOAD 



is dictated by the 



desired arc operating condition, improved electrical sta- 
bility at each condition can be effected only by increasing 
R^. This, in turn, necessitates an increase in the open- 
circuit voltage. Accordingly, for stable arc operation, 
it is generally advisable to use the highest possible 
settings of both the ballast resistor and the open-circuit 
voltage of the power supply. 



positive slope portions of the arc generator load lines . 
Nevertheless, the maximum open-circuit voltage capability 
of the power supply can place a restriction on high cur- 
rent operations. For example, if all of the ballast 
resistance has been removed in order to operate in condi- 
tion-a (Figure II. 3), further operation at higher arc 
currents (position-a ' ) or operation at greater values of 
ifi and/or L (position-a") can be achieved only by increasing 



* Note that operation at the initial point of intersection 
(e.g., points x and y in Figure II. 3) is not likely since, 
according to Equation II. 2, the arc would be inherently 
unstable at this condition. 



Irrespective of the value of R^, Equation II. 3 will 
be satisfied as long as operation is confined to the 



* 
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the open-circuit voltage. Therefore, aside from possible 
over heating of the electrodes, high current operation 
is limited by the maximum obtainable open-circuit voltage 
of the power supply. 

According to the stability requirement expressed by 
Equation II . 3 , operation on the negative slope portions of 
the arc generator load lines (e.g., position-b) is in- 
herently less stable than operation in the positive slope 
regions. In addition, a lower limit for stable operation 
exists at the current for which the slopes of the load and 
supply characteristics are equal (e.g., position-c) . In 
practice, the arc will extinguish prior to this lower 
limit; however, lower current operation (position-c') or 
operation at higher valves of ifi and/or L (position-c") can 
be achieved by increasing the open-circuit voltage and 
increasing R^. This procedure can be continued until 
the maximum open-circuit voltage of the power supply (or 
the maximum value of R^) is reached. 

For the duct lengths and operating gas flow rates 

employed in this study, the point at which the slope of 

the load characteristic is equal to zero ( ( - 5 -=-) = 0 ) 

dI LOAD 

occurs at an arc current of approximately 25 amp for Ar , 
150 amp for N^, and at a current in excess of 400 amp in 
He. Therefore, according to the previous arguments, low 
current operations in all three gases and high current 
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operations in N 2 and He are limited by the power supply. 
High current operation in Ar is limited by the maximum 
power that can be dissipated by the main ballast resistor. 
It is further noted that over heating of the anode (as 
evidenced by ablation of the anode surface) occurred at 
the higher arc currents in (I > 200 amp) and in He 
(I > 100 amp) . 

II. 3 Diagnostic Methods and Related Instrumentation 
The overall and local arc-column characteristics 
are determined by measuring the potential difference 
between the "floating" duct segments and the cathode. The 
arc voltage measurements are used both to determine the 
axial distribution of the electric field intensity and, 
along with the arc current and calorimetric measurements, 
to determine the mixed-mean enthalpy of the gas at each 
segment location. The arc voltages are measured with a 
high input impedance (4 megaohm) strip-chart recorder 
(Electronik-19 by Honeywell) . The arc current is deter- 
mined by monitoring the voltage drop across a calibrated 
shunt resistor with a digital voltmeter (Dana, model 
# 54030). As shown schematically in Figure II. 4, further 
circuitry is available to permit up to 400 ma of current 
to be drawn from a duct segment while its potential (with 
respect to the cathode) is monitored. 



29 





Figure II. 4 Circuit Used for Impressed-Current (I a ) 
Method of Wall Potential Measurement 
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The energy loss to the duct segments, to the 
cathode, and to the remaining arc generator components is 
determined calorimetrically . Accordingly, each of these 
components is equipped with a specially designed tempera- 
ture difference (AT)- transducer. The AT transducers 
(Figure II. 5) employ a single pair of copper-constantan 
thermocouples, mounted in plexiglass inlet and outlet 
plenum chambers, and wired in series to provide a voltage 
signal which is directly proportional to the coolant 
temperature rise. Considerable care was taken in the 
shielding and calibration (against mercury-in-glass 
thermometers) of these transducers and they are accurate 
to within 5% of the indicated reading. The transducer 
output voltages are amplified (Dana model #3850-V2) prior 
to being displayed on the digital voltmeter. The abso- 
lute temperature of the operating gas at both inlet and 
discharge locations is measured with a copper-constantan 
thermocouple (as referenced to a copper-constantan ice- 
junction) . The gas temperatures are monitored by the 
digital voltmeter. 

A unique radiative heat transfer gage is used to 
delineate the radiative contribution to the total wall 
heat flux. As shown in Figure II. 6, the radiation gage 
consists of two semicylindrical copper segments which 
are thermally isolated from one another, individually 
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Figure II. 5 AT-Transducer Assembly 
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Figure II. 6 Radiation Gage 
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water cooled, and supplied with individual AT-transducers . 
The inside wall of each half segment is coated with either 
a highly reflecting (aluminum) or highly absorbing (black) 
surface. The overall dimensions of the assembled gage 
are precisely the same as those of a standard duct seg- 
ment. A thorough discussion of the principle of opera- 
tion, the fabrication, and the associated experimental 
problems involved in the use of the radiation gage are 
left to the Appendix. 

The operating gas flow rate and the coolant flow- 
rates through the cathode and those duct segments serving 
as calorimeters are measured with precision accuracy 
(± 1%) flowmeters (Fischer & Porter model no. 10A0735) . 

The two gas flowmeters have a combined range of from 0.02 
to 1.50 lb/min in argon. At the time of each flowmeter 
reading, the inlet gas temperature and pressure are 
measured and used to correct the indicated flow rate 
should they differ from the temperature and pressure of 
calibration. An additional precision accuracy flowmeter 
is used to measure the total coolant flow rate to the 
plasma generator . 

With measurement of the arc current (I) and axial 
voltage distribution (E) , the gas inlet temperature ( t j NLET ^ 
and flow rate (m) , and the power losses to the duct 
segments (PLS) and to the cathode (PLC) , the axial 
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distribution of the gas enthalpy can be determined by 
performing energy balances at successive axial locations. 




Figure II. 7 Energy Balance Diagram 



With reference to Figure II. 7, the enthalpy of the gas 
at the exit plane of a segment in the arc heating region 
is given by 



j 

E . I - E (PLS) . - PLC 
3 i=l 1 



h . = C T 
3 p INLET 



+ 



A 



(II. 3) 
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This value is then algebraically averaged with the enthalpy 

determined by the previous ( j — 1 ) energy balance, and 

the resulting "mixed-mean" enthalpy is assigned to an 

t h 

axial location corresponding to the center of the j 
segment . 

To experimentally verify the accuracy of the diag- 
nostic procedures employed in this study, an attempt was 
made to "close off" the energy balance for the entire 
plasma generator. While the precision of this energy 
balance is limited by unavoidable ambient losses, it was 
still possible to obtain closure to within the cumulative 
experimental uncertainty of the component measurements 
(± 20%) . A more precise check is one based on a local 

energy balance performed on a segment located in the 

* 

fully-developed portion of the arc heating region . In 
every case, closure of this local energy balance was 
obtained to within 5% - a value comparable to the total 
uncertainty of these measurements (±5-10%) . The consistent 
reproducibility of data taken over a two year period pro- 
vides further confidence in the accuracy of the measure- 
ments as well as the consistency of the overall system. 

A detailed uncertainty analysis is used for estima- 
ting the uncertainties associated with each measured 



* The nature of this energy balance will be discussed 
in detail in the following chapter (Section III. 3). 
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variable and for calculating the propagation of these 
uncertainties into each reported result. The method used 
in this analysis is the standard approach for single-sample 
experiments . The uncertainty calculations, as well 
as the conversion and reduction of raw data, are per- 
formed on the computer. The thermo-physical properties 
required for determining the mean temperature, T m =T (h^) , 

and the Nusselt number-Ohmic heating parameter correlations 

[371 

(k and a ) are obtained from DeVoto and Drellishak, 

m m 

[ 3 g I r 331 

Knopp, and Cambel 1 in the case of Ar; Lick and Emmons 

and DeVoto and Li^^ in the case of He; and Drellishak, 

. [4H [421 

Aeschliman, and Cambel and Avco for N2* The re- 

ported uncertainties of these properties are also consi- 
dered in the aforementioned uncertainty analysis. A more 

detailed description of the arc facility, diagnostic tech- 

[431 

niques, and uncertainty analysis is given by Lukens 
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CHAPTER III EXPERIMENTAL RESULTS 

III.l Range of Measurements 

Experimental data were acquired for Ar , N^ , and 
He arcs over the following ranges of laminar flow operating 
conditions : 

Argon; I = 35 - 240 amp, ih = 0.03 - 0.11 lb/min, 
constrictor length = 5 in (20 segments) . 

Nitrogen; I = 100 - 200 amp, ifi = 0.03 lb/min, 

constrictor length = 5 in (20 segments) . 

Helium; I = 50 amp, m = 0.07 lb/min, 

constrictor length = 2.5 in (10 segments) . 

All measurements were made with a 1.0 diameter copper 
duct operating at near atmospheric pressure. The assertion 
of laminar flow is based upon the flow transition studies 
by Incropera and Leppert L ^ J and by Runstadler in 

which the onset of turbulence is reported to occur at a 
mean Reynolds number (based on properties evaluated at T ) 
of 700 to 900. The largest value of the mean Reynolds 
number encountered in this study is approximately 600 
(Argon, I = 35 amp, m = 0.11 lb/min). Hence, the above 
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stated ranges of operating conditions are assumed to be 
representative of laminar flow. 

Rather than presenting a complete tabulation of the 
data acquired for the inlet region, certain data were 
selected to demonstrate the effect of varying gas flow 
rate and arc current on the development characteristics 
of the arc. All of the acquired data pertaining to the 
asymptotic region are presented along with related data 
from other investigators. The interpretation of the 
experimental results as well as comparisons with related 
theories will be the subject of Chapter IV. 

III. 2 Electrical Characteristics 

Prior to presenting the experimental data for the 
arc-column (electrical) characteristics, a discussion is 
in order concerning the method employed in the acquisition 
of this data. 

A method frequently used in the measurement of 
arc-column characteristics involves use of the duct seg- 
ments as "floating" potential probes . It 

is commonly recognized, however, that a floating duct wall 
will assume a potential which is negative with respect 
to the plasma core by an amount equivalent to the plasma 

i k 

sheath potential . Accordingly, on the assumption that 



* Further treatment on the formulation and characteristics 
of plasma sheaths is given in references [30] and [45] . 
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the sheath potential is appreciable, some investigators 
have questioned the suitability of using a floating duct 
wall to measure plasma potentials. 

In reference to their cascade arc potential measure- 
ments, this question has been raised most strongly by 

r 6 7i 

Pfender, et al . ' . In performing their measurements, 

they allowed a current (I & ) to be drawn from the duct 
wall while they simultaneously monitored the potential of 
this wall with respect to the cathode. These measure- 
ments indicated that the recorded wall potential was 
strongly influenced by the magnitude of the current drawn. 
In particular, they observed that the voltage increases 
by 40 to 100% over the floating wall (I = 0) potential 
when a current of up to 50 ma is drawn. Above 50 ma, the 
segment potentials were observed to be independent of the 
current drawn. With regard to this observation, they 
argued that, for I > 50 ma, the charge separation which 
otherwise would exists in the plasma sheath region is 
exactly cancelled and, consequently, the duct wall is at 
the potential of the plasma column. They therefore con- 
cluded that floating wall potential measurements do not 
provide a suitable indication of the actual plasma 
potential, but that accurate measurements of this para- 
meter can be made by recording the wall potential while 
a current in excess of 50 ma is drawn from the wall. 
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In contrast, Emmons indicates that the sheath 
potentials for an argon-copper duct combination are 
approximately 1 to 2 volts, for which case the floating 
wall potentials are suitably close to those in the arc 
column. In support of this contention is the experimental 

[15] 

technique used by Maecker wherein every attempt was made 

to eliminate impressed wall currents (I ) while performing 

3 . 

* T461 

plasma potential measurements . Schreiber , on the 
other hand, has measured significant voltage increases 
with current drawn, but he attributes these increases to 
voltage drops across the relatively cool gas boundary lay- 
ers. Additional support for Schreiber 1 s deduction can 

be inferred from the measurements made by researchers at 
[32 331 

Dartmouth ' . While the Dartmouth data were used 

primarily in support of anode studies, the arc configura- 
tion and external circuitry used for these measurements 

are similar to those employed in this study (Figure II. 4) 

[71 

and to those used by Pfender, et al. . 

In an attempt to resolve these inconsistencies, the 
variation in the duct wall potential with current drawn 
was measured at various axial locations and arc currents. 

As shown in Figure III.l, these measurements indicate that 



* Maecker contends that the measured wall potential 

difference between two segments will be in error by an 
amount proportional to the impressed current, I . 

cl 
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Figure III.l Dependence of Wall Potential on Curren 
Drawn, Arc Current, and Axial Location 
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the percentage increase in wall potential with increasing 
I diminishes both with increasing axial distance and 
with increasing arc current. They also indicate that 
the wall potential increases continuously with increasing 
I over the entire range of I values investigated. In 

cl 3 . 

this respect, the distributions of Figure III.l are 

qualitatively similar to the potential distributions 

[15] 

recorded by Maecker and to the lower current ranges 

of anode conduction reported by Runstadler . 

Simultaneous calorimetric and wall potential mea- 
surements indicate that the total wall heat flux to a 
duct segment does not vary appreciably (< 3%) over the 
entire range of I . In addition, the value of the 

cL 

asymptotic electric field intensity (E ) was observed 

X / 3 

to be independent of the current drawn. This fact 

[11] [32] 

is discussed by Emmons and Land and by McKee, et al . 

and is strongly supported by the experimental measure- 

[141 

ments of Morris, et al. , wherein the voltage gradient 
was measured using three independent methods . With the 
first method , the value of E was computed from the 

X , 3 

asymptotic energy balance equation (this equation is dis- 
cussed later in this chapter) along with calorimetric 
measurements of the power losses to the duct segments. 

The second method, which is similar to that employed by 
Maecker, utilizes two of the duct segments as potential 
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probes. Unlike Maecker ' s method, however, a small current 
was supplied to the arc through these probes (in an attempt 
to cancel the voltage drop across the "contact resistance" 
between the probes and the arc column) while the potential 
difference between the two duct segments was measured. 

This potential drop was then divided by the distance be- 
tween the two probe segments to give the value of . The 

third method used by Morris, et al . involved a cathode 
withdrawal experiment wherein the total arc current was 
held constant while the total (cathode to anode) voltage 
drop was recorded as a function of the total (cathode 
to anode) duct length. These three methods, when applied 
to the central region of a hydrogen arc, resulted in values 
of the electric field intensity which agreed to within 
approximately 10% - a value comparable to the combined 
uncertainties of these measurements. 

The agreement demonstrated by these experiments is 
presumably a result of the conditions which characterize 
the fully-developed region. That is, regardless of the 
magnitude of possible radial electric field components, 
the voltage difference between any two axial points in 
the fully-developed region is invariant as long as the two 
points are at the same radial location. If this were not 
the case, the axial component of the electric field 
intensity would vary with axial location and the flow 
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would no longer be asymptotic. The independence of E 

X f cl 

with I was tested under the most severe condition 

cl 

(I = 50 amps, I = 400 ma) of this investigation and, as 
shown in Figure III. 2, the slopes of the linear portions 
of the wall potential distributions (E ) are the same 

X f cl 

for I =0 and I = 400 ma. It is further noted that the 

cl Si 

electric field intensity distribution on either side of a 
segment from which current is being drawn is unaffected 
by the wall potential variations of that segment. 

Although it is concluded that the current drawn 
from a duct segment (up to I = 400 ma) does not apprecia- 

Si 

bly effect the total wall heat flux or E , the increase 

X f Si 

in wall potential with increasing I will nonetheless have 

a direct effect on the energy balance measurements used 

to determine the mean enthalpy of the gas, particularly 

at the lower arc currents. As indicated in Figure III. 3, 

the mean enthalpy of a 50 amp, argon plasma is considerably 

greater when computed using wall potentials measured with 

I a = 400 ma than it is when computed with floating wall 

potentials. While the relative accuracy of these two 

distributions could best be determined by comparison with 

experimental enthalpy profile data, such data are not 

available. However, a comparison is made with the non- 

1251 

equilibrium numerical solution by Clark . This 
solution, as applied to the asymptotic region of an 
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Figure III. 2 Axial Distribution of Wall Potential with Current Drawn 
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Figure III. 3 Axial Distribution of Mixed-Mean Enthalpy with Current 
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atmospheric argon plasma in a 1.0 cm diameter duct, pre- 
dicts electric field intensity and total wall heat flux 

values which agree with measured values to within the 

* 

experimental uncertainties (< 6%) . As indicated in 

Figure III. 3, this solution also predicts an asymptotic 
mean enthalpy which agrees well with the experimental 
value determined from the use of floating wall potentials 
While the conclusions that can be drawn from this compari 
son are, at best, tenuous, the agreement between theory 
and experiment shown in Figure III. 3 indicates that the 
floating duct wall serves as the more suitable potential 
probe . 

After modifying the wall boundary conditions to 

include current drawn from the duct wall, the nonequili- 

F251 

brium solution by Clark (which includes a plasma 

sheath model) was used to investigate the effect of I 
on the plasma sheath potential. At I = 50 amp and 
I = 400 ma , the nonequilibrium solution predicts a 
decrease (from its floating wall value) in the plasma 
sheath potential from 3.60 volts to 1.75 volts. This 
predicted reduction in sheath potential is considerably 
less than the increases in wall potential measured under 



* 



Comparisons between experiment and the nonequilibrium 
model are presented in Chapter IV. 



